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Abstract 

m : 

. We prove a uniqueness theorem for an inverse boundary value problem for the 

Maxwell system with boundary data assumed known only in part of the bound- 
On I ary. We assume that the inaccessible part of the boundary is either part of a plane, 
■ or part of a sphere. This work generalizes the results obtained by Isakov [I] for the 
. Schrodinger equation to Maxwell equations. 

Introduction. Let f2 C be a bounded domain witli C^'^ boundary, and let e, /i, a be 
functions in f2 [e is the permittivity, /i the permeability, and a the conductivity). We 
will assume that the coefficients satisfy the positivity conditions 

> 

(N 
O 

Tj- ■ Let D = —iV, let u be the exterior unit normal to dQ, and consider the time-harmonic 

CN ■ Maxwell equations for the electric field E and magnetic field H in Q, 
O 

g: lDAH + u;^E = 0, 

\ DAE-LjfiH = 0, > 

^ ' with the boundary condition 

c3 ■ u A H = a on dVt. (0.3) 

This is the magnetic boundary value problem for the Maxwell equations. Here we use 'A' to 
denote the vector product in M^, and V A F is the curl of the vector field F. When posed 
in correct function spaces this problem admits a unique solution [E, H) when the angular 



7 = £ + ia/uj, e > 0, /i > 0, (J > inil. (0.1) 
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frequency lo > does not belong to a discrete set of resonant frequencies. The impedance 
map A(7,/i) is then defined by 

A(7,/i) : a^u A E\dn- 

In analogue with the Calderon problem of electrostatics (see [C], [N] and |SUj ). we consider 
the inverse problem of determining the electromagnetic parameters 7 and fi from the 
boundary measurement A(7,/i). Especially, in this work we assume that we can measure 
the values of A(7,/i)(a) only on a nonempty open subset F of dfl, and only for tangential 
boundary fields a supported in F. We further assume that the inaccessible part of the 
boundary Fq = dQ \ F is either part of a sphere, or part of a plane. Similar results were 
proved by Isakov ^ for the inverse boundary value problem for the Schrodinger operator. 
Our work extends this method to Maxwell equations. 

Let us now formulate the main results of this paper in a precise way. Given Q as above, 
we define the spaces 

Hl,^{n) = {Fe H\nf; Div(z/ A F) e H^/^dQ)}, 
TH^^.^idn) = {ae H^'\dnf- z/ ■ a = 0, Div(a) G H^'\dn)}, 

where Div denotes the surface divergence on dVL (see for example [CP] and [_OPSlj for more 
information). Provided that /i and 7 are coefficients in C^{VL) which satisfy (10.11) . then for 
uj outside a discrete set of resonant frequencies, the magnetic boundary value problem 
(irr^ has a unique solution in H}^-„{n) x H^^^^iVt) for any a G TH^{l{dn) (see |UPSl] l 
The impedance map is well defined, and 

A{^,f,):TH'^{l{dQ)^TH'£{dQ). 

Our first main result is the following. 

Theorem 1. Assume that f2 c {M^; < 0} is a C^'^ domain. Let Tq = d^l D 
{x3 = 0} and F = dQ \ Tq. Assume we have two pairs of electromagnetic parameters 
iljyf^j) ^ C"'(fi) X C^{Q), j = 1, 2, satisfying the positivity conditions (10. ip and the 
following boundary conditions: 

7i = 72, /^i = /^2 up to order one on T, and (0.4) 

{there exist extensions of 'jj and nj into M.^ which are , , 

invariant under reflection across the plane {x^ = 0}. 

Suppose that uj is not a resonant frequency for (71, /ii) or (72, /i2)- If the impedance maps 
restricted to F coincide, 

1 /2 

A(7i,/xi)(a)|r = A(72, /i2)(a)|r for all a G THj^^^{dQ) with supp{a) C F, 
then the electromagnetic parameters are equal, i.e. 71 = 72 and fii = ^2 in VL. 
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Before formulating the second main theorem of our paper, we want to comment on the 
assumptions in Theorem 1. The condition flU.4p is not important, and we expect that a 
suitable boundary determination result would allow to remove this condition completely 
(in the case of C°° coefficients such a result is given in | JM] ) . The condition (10.51) comes 
from the method of proof, which requires that the Maxwell equations are invariant under 
reflection across {xs = 0}. The following is a necessary condition for (lU.Sp to hold: 

For the sake of definiteness, we have stated Theorem 1 in terms of the impedance map. 
However, the proof of Theorem 1 extends to the case of restricted Cauchy data sets 

Cr{-f, /i) = {(z/ A E\r, v A E\y)\ {E, H) G H'^i^f x H'^i^f solves 

In this setup Vt can be a Lipschitz domain and > any number, and the assumption 

Cr(7i,/Wi) = C'r(72,/i2) 

for two pairs (7j,/ij) satisfying the conditions in Theorem 1 will imply that 71 = 72 
and /ii = ^2 in ^- This result allows domains with transversal intersections of dVl and 
{xa = 0}, in which case (10. 6 p is also sufficient for (10. 5p . It is well known (see for example 
|MM] , |Slj and |S2j and the references given in these articles) that Maxwell equations may 
not admit if ^-solutions even with boundary data in H^^'^^dVl) if the domain is non convex: 
the solution may have a finite dimensional singular part near the conical singularity, so the 

1 /2 

image of the impedance map might be slightly larger than just THjj-^^{dfl). However, all 
solutions needed in the proof are restrictions to Q of if^-solutions defined in a neighborhood 
of Q (actually they are defined in the whole M^), and hence it is enough to restrict only to 
the part of the graph of A(7,/i) that is actually needed. 

Let us now formulate our second main theorem: 

Theorem 2. Assume Q C Bq is a C^'^ domain included in an open ball Bq of M.^ of 
positive radius, and let Tq = dQ fl OBq, Tq 7^ OBq, T = dQ \ Tq. Assume we have two 
pairs of electromagnetic parameters {'~fj,fij) G C'^(fi) x C^(fi), j = 1, 2, satisfying the 
conditions of Theorem 1 hut with respect to OBq instead of {x^ = 0}. If uj is not a resonant 
frequency and if the impedance maps restricted to T coincide, 

1 /2 

A(7i, /ii)(a)|r = A(72, yU2)(«)|r for all a G TH^-^^{dVt) with supp{a) C F, 

then the electromagnetic parameters are equal, i.e. 71 = 72 and fii = ^2 ^■ 

The unique recovery of electromagnetic parameters from the scattering amplitude was 
first proven in |CP] under the assumption that the magnetic permeability /i is a constant. 
The unique recovery of general C^-parameters 7 and /i from full boundary data was then 
proved in [OPSlj , and later simplified in |OSj , and the uniqueness for the inverse scattering 
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problem was proved in |Saj . Boundary determination results were given in |McDlj and 
[JM] . and the more general chiral media was studied in |McD2j . For a slightly more general 
approach and more background information, see also the review article |0PS2] . Note 
that in this work we need, for technical reasons, two more derivatives of the parameters 
compared to the full data problems. 

Inverse problems with partial data for scalar elliptic equations have attracted consider- 
able attention recently. There are two main approaches: the first uses Carleman estimates 
|BU] . |KSU] and the second uses reflection arguments [I]. In the first approach Q can be 
any domain but one needs to measure part of Cauchy data in a small set F C dQ and the 
other part in a neighborhood of dQ \ F. In the second approach it is enough to measure 
Cauchy data on a fixed small set F, but one has the restriction that dQ \ F has to be part 
of a plane or sphere. 

In this paper, we apply Isakov's reflection method [I] to the Maxwell system to prove 
Theorems 1 and 2. As far as we know, these are the first partial data results for an elliptic 
system which are analogous to [I] (an analog of [KSUj for a Dirac system was recently 
proved in |STj ). The strategy of the proof is similar to [OPSl] and [OSJ: we construct 
special solutions to the Maxwell system, insert them in a suitable integral identity, and 
recover the coefficients from an asymptotic limit. 

However, several issues arise when trying to combine this strategy with the reflection 
method. The integral identities in |UPSlj and |OSj do not seem to go well together with 
partial data, and the final identity which we use for recovering the coefficients is new. 
In this identity, many terms survive in the asymptotic limit, and one needs to carefully 
manipulate these terms to determine the coefficients. Also, there are cross terms involving 
products of the original and reflected solutions. In [I] these cross terms were handled by 
the Riemann-Lebesgue lemma. In our case we need to prove a substitute for the Riemann- 
Lebesgue lemma, and here the C'^-regularity of the coefficients is used. 

The structure of the paper is as follows. In the first section we augment the Maxwell 
equations into a Dirac type elliptic system and prove a crucial factorization result for 
a second order matrix Schrodinger operator. This follows in principle |USj . but with 
somewhat different notations and some simplifications. In order to keep the article self- 
contained, we have included complete proofs. In the second section we construct the 
complex geometric optics solutions (CGO solutions for short) needed in the proof. In this 
section, following Isakov [I], we also perform the reflection across {0:3 = 0} and analyze 
the behavior of CGO solutions for Maxwell equations under this operation. In the third 
section we prove the integration by parts formula needed in our argument, and in the 
fourth section we perform the asymptotical computations needed in our proof. Theorem 1 
is then proved in section five, and finally using the Kelvin transform, the proof of Theorem 
2 is reduced to Theorem 1 in the last section. Note that with the exception of the last 
section, we use the vector field formulation for our equations. Only when analyzing the 
behavior of Maxwell equations under the Kelvin transform does it become useful to revert 
to a formulation using differential forms. However, we feel that the major part of this work 
is easier to read when written in vector field notation. 
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1. Augmented Maxwell system and reduction to Schrodinger equation. In 

this section we recall an argument used in [OS] . Roughly, the idea is to augment Maxwell 
equations in such a way that, with a suitable rescaling of the fields, the system can be 
transformed into a Schrodinger equation. This argument can be written in a slightly more 
general setting than the one described above, see |0PS2j . Let U C be an an open set 
along this section. Consider the time-harmonic Maxwell system in U given by 

( DAH + u^E = , . 

\ -D AE + ujfxH = ^ ^ 

where > is a fixed angular frequency, 'j = e + ia/u and fi,e,a E C'^iU) satisfy 

/i, e > 0, (T > in f/. 

Taking the divergence of both equations in (11. 7p we obtain the additional equations 

r D ■ i^E) = 
\ D ■ (fiH) = 0. 

These will be useful in writing the Maxwell equations as an elliptic system. Namely, 
denoting a = log 7 and /3 = log fi we may combine these into four equations 

D ■ E + Da- E = 

-DAE + uj^H = 
D ■ H + DP ■ H = 
D AH + LU-fE = 0. 

This gives rise to the 8x8 system 
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= 0. 



Here * means that in these positions we may insert anything and would still get the same 
equations. 

We wish to consider 8- vectors X* = ( $ if * | \& -E* ) , where and $ are additional 
scalar fields. We also want the system to be elliptic, in fact, we want the principal part to 
be a Dirac-type operator. To this end, consider two 4x4 operators (acting on 4-vectors 
( $ if* ) for instance) 
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These operators satisfy 



P+{D)P_{D) = P_{D)P^{D) = -Ah, 
P+{D)* = P-{D), P_{Dy = P+{D). 



Motivated by the Pauh-Dirac operator in R^, we choose the principal part to have the 
block form 

P_{D) 
P+iD) 



Observe that P^{D)P^{D) = -A/g and P^{Dy = P^{D). The following notations will 
also be used in the rest of the work. Let us denote 



P^{a,b) = 
where a, 6 e C^; and 






P^b) 


P+(a) 






P^a,b) 








P_(a) 






(1.8) 



diag(^,S) = 



A 








B 



where A and B are 4x4 diagonal C-matrices. When A — A1/4 and B — A2/4 with 
Ai, A2 G C let us write 

diag(Ai, A2) = 



A1/4 








A2/4 



Note that P^(a, 6)P^(6, a) = diag(6 • 6, a • a) and P^{a,b)P^{b,a) — dia,g{b ■ b, a ■ a) for any 
a, 6 e C^. We shall write for short P^(a) and P^(a) when b — a, hence P^{a)P^{a) — 
{a ■ a)Is and P^{a)P^{a) = {a ■ a)Is for any a G C^. Prom these last equalities follow the 
anticommutation formulas 



P^(a)P^(6) = -P^(b)P^(a) + 2(a • 6)78 
P^{a)P^{b) = -P^{b)P^{a) + 2(a • 6)78, 

for a, 6 e C^. Moreover, the commutation formula 

P^{a)P^{b) = P^{b)P^{a) 



(1.9) 
(1.10) 



(1.11) 



holds for a,b E C^. This formula follows easily from P+(a)P+(6) = (P+(a)P+(6))* 
P-{b)P-{a). More generally 



P^{a,b)P^{c) = P^{c)P^{b,a), 



(1.12) 



for a,b,ce C^. 

Por the potential there is seemingly considerable freedom in the choice of the entries 
marked with *. However, we later wish to reduce the Maxwell equations into a matrix 
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Schrodinger equation with no first order term so a proper choice of the extra entries will 
be crucial. We will take 



V, 



/ u}^ Da- \ 

couls Da 

D(3- u-f 

\ Dp u-fh ) 

The augmented Maxwell system is then 

(P^(D) + V;,,^)X = in[/. 

The solutions of this system for which $ = = correspond to the solutions of the 
original Maxwell system in f/. 

We may write the potential V^_^ above in the form 



V^,^ = a;diag(/x, 7) + 1{P^{D(3, Da) + P^{D(3, Da)). 
Next we rescale the X as follows. Let 

be defined by 



:i.i3) 



Then 
where 



Y ^dieig{ii'/\j'/^)X. 
(P^D) + V,,,)X = diag(7-^/^ /.-V2)(PT(^) + w,,,)Y 





\P+{Da) 




KI4 



^Kls + -P^{D(3,Da) 



and K = a;(7/i)^/^. Here we used that 



P^(Z^)(diag(/x-l/^7-'/')>^) = diag(7-^/^/x-^/^)[P^(D) - \p^{D(3, Da)]Y 

Hence 

{P^{D) + V^,^)X^O ^ {P^{D) + W^,^)Y^O 
when X and Y are related as above. This scaling is motivated by the following result. 



Lemma 1.1. One has 

(P^iD) + W,,,)iP^iD) - Wl^) = -AIs + Q,, 
{P^{D) - Wl^){P^{D) + W,,,) = -AIs + Q',, 
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where the matrix potentials are given by 

/ Aa 

2V'^a-Aal3 











2Dk- 
\ 2Dk 



A/3 
2V2/3- A/373 y 
2Dk- \ 
2Dk 



and 



1 


/A/3 
2V2/3- A/3/3 





2 


[ 


Aa 




2V'^a-Aal3 






2DkA 




-2DkA 





I 

^lth^''j = {dl^^,J)l^^,. 

Proof: We have 

{P^{D) + W,,,){P^{D) - Wl^) = -Ah + W,^,P^{D) - P^{D)Wl^ - W,,,W\ 
We claim that 



t 



W,^,P^{D) - P^{D)Wl^ 



Aa 
2V2q; - Aq;73 








A/3 





2V2/3- A/373 



\ 



This can be proved by a direct calculation. However, we will use symbol calculus which 
gives a slightly more elegant proof. In standard (left) quantization, since the symbol of 
P^{D) is P^(0 and since d^.P^i^) = d^.W^,^ = for j = 1, 2, 3, we have 

3 

w,,,p^D) - p^{D)wi^ = op{w,,,p^{o - p^mi, - J2 duP^m^^wi,) 

m=l 

3 

m=l 
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The second equality holds because 



W,^,P^{i) - P^iml, = \p^{D^,Da)P^{i) - ^P^iOPHDa,D(3) = 0, 
which is true by (11.121) . Consequently 



W^^^P^iD) - P^{D)W. 



op(p^(DX^) 




1 f P_{D)P_{Da) 





2 V 


P^{D)P^{D(3) 



On the other hand 







P^{Dk) 





-P^{Dk) 



since 



-{Da + Dp) = Dlogifi-fY^^ = k-^Dk. 



This shows the desired form of Q^^-^. 

It remains to prove that Q'^ ^ has the expression given. We have 

{P^{D) - Wl^){P^{D) + W,,,) = -AIs - Wl^P^{D) + P^{D)W,,, - Wl^W,,,. 

Note that W!, ^, = W^^^. Hence, changing the role of and 7 we obtain 



-Wl^P^{D) + P^{D)W^,^ = -W^,^P^{D) + P^{D)W'^ 



P^{D)P.{DP) 








P+{D)P+{Da) 



+ P^{Dk). 



and 



{^^ + 1{DPY)I, 


P+{Dk) 


P-{Dk) 


{K' + i{Day)h 



This shows the desired form of Q' ^. 



□ 



The crucial point in the above lemma is of course the vanishing of the first order terms. 
Also, it is worth pointing out that we will need both second order operators given in Lemma 
ll.ll in our work. 

2. Construction of CGO solutions. In this section we construct CGO solutions for 
the rescaled Maxwell system and then we produce solutions for the original Maxwell sys- 
tem. In order to get such solutions we apply the Sylvester-Uhlmann method (given in [SUj) 
to the matrix Schrodinger equation obtained in Lemma fTTTl This argument was introduced 
in [OS]. 
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For the sake of completeness we recall the Sylvester-Uhlmann estimate whose proof can 
be found in (SU], §\- Let g G L°°(M3; C) with compact support and -1 < 5 < 0. If C G 
such that C ■ C = k'^ with \(\ large enough, then for any / G there exists a unique 

u G H] solving 

{A + 2C-D + q)u = f. 

Moreover 

MHi<c\cr'\\fh.^^, 

for < s < 2 with C a constant independent of ( and /. The norms are given by 

\\f\\%= [ {i+\xmf{x)\'dx 

and 

= ||(1 + 1^1 ) ^ ^||//»(R3) ■ 



Assume that we can extend and 7 to the whole space so that for some positive 
constants and /iq we have j — Eq, fi — fio & Cq(]R^), and let Q^^^ = u'^EqUqI^ + Q^^^ which 
is thus compactly supported. Denoting k = c<j(eo/Uo)^^^ we are thus led to the equation 

{-{A + e)Is + Q,,,)Z = (2.14) 

Let ^ G M'^ any non-zero vector. Consider G such that ( ■ ( = k"^ and 

where r > 1 is a free parameter controlling the size of (, and C = ''7i + ^^2 with r]i,ri2 G M^, 
Vi ' V2 = and \rji\ = \rj2\ = 1. Let Zq = Zo{() be a vector which does not depend on x 
and which is bounded with respect to the parameter r. 

Proposition 2.1. Assume that -1 < 6 < 0, and considers > such that -1 < 6+6 < 0. 
There exists a CGO solution of the form 

Z = e'^-\Z^ + Z_i + Z,) 

for the equation \2.1J^ in M^, such that 



for < |a| < 2. Moreover, if we denote 

R = lim rZ_i M = lim Zq, 

r— >+oo r— >+oo 

then 

2{C-D)hR = -Q^^^M (2.15) 

in M?. 
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Proof: Let ( & he as above, note that 

e-<- (_(A + e)Is + Q,,^) e^'^^Zo + Z_^ + Z,) = 
= Qm,7^o + ((-A + 2C • D)Is + g^,^)(^-i + Zr) 

in M^. Let x be a cut-off function which is identically 1 on a neighborhood of {x G : 
\x\ < p} where p is so large that Q C {|x| < p}, and define 



(2vr) 



For r > let Xt{x) = xi^/'^) 



Z_, = —XriC-D)-\-Q,,,Zo). 
It 

It is well known that (C ■ D)^^ : is bounded with norm of order \C\~^ ( |SU] ) . 

hence 



■5 r 



"5+1 



Moreover, differentiating we get 



for all 1 < |a;| < 2. Here we need the C^-regularity of the coefficients. On the other hand, 
let Z^ be the solution of 

-((-A + 2C-/^)/8 + Q;.,^)^r = ((-A + e(r)-Z})/8 + Q^,^)Z_i + (C-/^Xr)(C-^)"'(-QM,7^o), 
satisfying 

WZrW^, < c|c|i"i-^ ((-A + e{T) ■ D)h + Q,,,)Z-, + (C • Dxr){C- D)-\-Q,,,Zo 



^S+1 



for < lal < 2. We then have the estimates 



4+1 



(X-Dxr)i(-D)-\-Q,,,Zo 



< Cr 



-2e 



{(■D)-\-Q,,,Zo 



o+e 



and 



-A + d{T) ■ D)hZ^ 



c 

< — 



S+1 



(C-D)-\-Q,,,Zo, 



2 


c 


H 





D{C-D)-\-Q,,,Zo) 



-A + e{r)-D){C-D)-\-Q^,,Z,) 



o+e 
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Let us explain the second inequality in more detail: consider e > such that — 1 < 6+6 < 0. 
Then 



iC-Dxr)iC-D)-'i-Q,,,Zo, 



[1 + \x\'Y^'\{C- Dxr){C- Dr\-Q,,,Zo)\'dx 



xgsupp(xT) 



< sup {{1 + \xn-^C ■ Dxrl''} I il + \x\'Y+^\iC-D)-\-Q,,,Zordx 



(C-D)-\-Q^,,Zo 



Therefore 



ID'^Z.^ 



r\\r? 



0{t 



for < |a| < 2. Finally, R satisfies the equation (12.151) by the definition of Z_i and the 
continuity of (C ■ D)~^ . □ 

Considering Z as in Proposition 12. II we have, by Lemma [1. 11 that 

Y = {P^{D)-Wl^)Z 
is a solution of the rescaled system 

(p^(D) + i^^,^)y = 0. 

Note that this solution can be written as 

Y = e'^^-^iYi + Yo + Yr) 

where 

Fi = P^(C)^o, 
Yo = P^OZ-i - Wl^Z^, 

Yr = (P^iD) - + {P^iD + C) - Wl^^)Zr, 



(2.16) 



(2.17) 



and 



Oir), \\D-Yo\y^^^ = Oil), WD'^YAl^.^u) = ^^(^''''"' 



(u) 



(2.18) 



for < |a| < 1 and U any bounded open subset of M.^. Recall that our aim is to construct 
CGO solutions for the original Maxwell system. As we mentioned above, Y will produce 
a solution of the original Maxwell system if = Y^ = 0. In order to get this condition 
we choose Zq in a special way. 



Lemma 2.2. // 

then we have 

for T large enough. 



((P^(C) - k)Zof = ((P^(C) - k)Zof = 

Y'={o {Y^y 1 (Y^y ) 
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Proof: This proof was given in [US] . Here it is included for the sake of completeness. 
As we mentioned above Y is solution of 

so it is also solution of 

(P^iD) - Wl^)iP^iD) + W,,,)Y = 0. 
From Lemma [1.11 we know that and are solutions of 



-(A + + g^y* = 

-(A + k^)Y^ + q^Y"" = 0, 



(2.19) 



where the potentials 



1 



g/3 = -\^(3 -K' + e- \{D(3)\ = -^Aa - + 1,' - -^{Da) 

are compactly supported. 

Adding i/c/s one gets ^ ^ 

Y = e'^-"(ri + Yo + Yr) 

where Yi = (P^(C) - kIs)Zo and % = P^(C)^-i + {klg - Wl^)Zo. This shows that 



(-A + 2(-D + + lo + Yrf = 0, 

(-A + 2C-D + + Yo + Yrf = 0. 



(2.20) 



It is known that these equations have unique solutions in L'j for |^| large enough (see |SU] ) . 
Since neither Y^ nor Y^ belong to unless Y^ = Y^ = 0, we choose Zq such that 

= ((^^(C) - kIs)Zof = 0, = ((P^(C) - ^18)^0)" = 0. 

In this way we have (Yi + Fq + ^r)*, (Yi + Fq + ^r)* £ Lj. Then, the unique solvability for 
\(\ large implies that F* = F"^ = for \(\ large enough. □ 



An explicit choice of Zq verifying the condition of last lemma may be done as follows. 



Let 



fC-a\ 
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\ ka J 




v* 
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kh 




\a) 




VC 


kh 


1 


\a) 



1 

r 



\{kh + \{PHO + P^{C))]m, 



(2.21) 



where we denoted m* = (O 6*|0 a* ) for any a, 6 G C^. The same choice was done in [OS] . 
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Next we perform the reflection argument given by Isakov in [I]. As in the statement of 
Theorem 1, let Q C {x^ < 0} be a C^'^ domain. With respect to the standard Cartesian 
basis e = {ei, 62, 63}, we introduce the reflection 

X = {xi,X2, X3) {xi,X2, -X3) =: x{x) 

and the reflected domain Q = {x{x) G M'^ : x G Q}. Consider 

o = nu int(ro) u n. 

By the assumption (10.51) . we can extend the coefficients fij and •jj into as functions 
which are even with respect to 0:3 for j = 1,2. We still denote the extended coefficients 
by II j and 7^. Note that these functions have the proper regularity to construct the CGO 
solutions as above. 

For the pairs (yUi,7i) and (/i2,72), for a given vector ^ = (^',^3) G satisfying 
l^'l > 0, we construct CGO solutions Zi,Yi,Xi and Z2, 1^25X2 in with complex vectors 
Ci and (2, respectively. We choose 

Ci = + ^ (r^ + ^) VI + {r' + k^) V2 (2.22) 

C2 = -\^-^ [r^ + ^) ^1 + {r' + k') (2.23) 

with r > 1 a free parameter controlling the size of | Ci | and | C2 1 , where rii,ri2 G verify 
I'?!! = 1^2! = I5 ^71 ■ ''?2 = and r/j ■ ^ = for j = 1,2. More precisely, for = (^1,^2,^3) = 
(f , 6) we choose 

1 ^1 
Vi = ^(^2, -^1, 0), ?/2 = ?/i A — = ^^(^^1^3, -66, l^'H- 

Note (j ■ (j = k"^ for j = 1,2. Let us denote 

( = lim — = lim — = ri2 + irji 

r— >+oo 7" 7" 

( = lim — = ri2 — irji 

and 

6(r)=2(Ci-rC) e2(r)=2(C2-rC). 
Observe that C,i{t) = C, + 0{t^^) and C,2{t) = — C + 0(r^-'^). Furthermore, if we denote 

Ri = lim rZi^ ^2 = lim tZ'^^ 

r— >+oo r— ++00 

Ml = lim M2 = lim Z^ 

r— >+oo T— ++00 
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then the equations 



2(C • D)IsRi = -g^i.^iMi 



hold in 



(2.24) 
(2.25) 



We now fix another orthonormal basis / = {/i, /2, fs} with 

/2 = ^(6,6,0), /3= (0,0, 1)^63, /1-/2A/3. 

With respect to the basis / one has 

e = (0,1^1,6)/, 
r^i = (l,0,0)^ = /i, 
V2^\^\-\0,-^3,W\)f- 

Therefore, 

G = + ^) f - (T= + *')"'|. I + (-^ + 

It is obvious that for any x,y gM.^ 



(2.26) 



X 



y = ^^jyj = ^^jy. 



f 



j=i i=i 

where x = {xi,X2,X3) = (x{,x|,x{)/ and y = (t/i, ^2, Z/3) = (z/i , 2/2> Z/a)/ with X3 = x{ and 

We next show that the reflection argument allows to construct CGO solutions whose 
tangential magnetic fields u A H on dQ have support contained in F. To this end we 
introduce 

Xj{x) = diag(-/4,/4)^j(i(x)) 



for j' = 1,2, where 



/I \ 
10 
10 

\0 -1/ 
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It is an easy matter to check that Xi and X2 are solutions of 

in Q (but also in Cl). Therefore, for r large enough, 

/ \ / \ 







Xi = Xi + X 



1) 







X2 + X2 



are solutions of the same equations and z/ A i^jlro = for j = 1, 2. Let us state the result 
obtained as the following proposition. 

Proposition 2.3. Corresponding to pairs (/ii,7i) and {^21^2)' '^'^'^ /^^ ^ vector 
^ G satisfying |^'| > 0, i/iere are CGO solutions Xi and X2 of the augmented Maxwell 
equations in VL such that 



iP^iD) + V,,,^^)X2 








where the terms X,- satisfy 



Xi = diag{^^ ^^'^ ^^'^)^\ + diag[-I^,l^Yi{x)) 



X, 



(i^a^(/i2^^^72^^^)(^2 + diag{-ii,i^Y2{x)) 



and the rescaled fields Yj are given by (12.161) and (I2.17P for complex vectors Q as in (I2.22p 
and (I2.23p . with ^ = Ci — (2- Moreover, for t large enough the fields Xj will solve Maxwell 
equations in Q and the tangential components of their magnetic fields will vanish on Tq. 



3. The orthogonality identity. We assume that uj is not a resonant frequency, so the 
Maxwell equations (10.21) have a unique solution {E, H) for any prescribed value of u A H 
on dfl as discussed in the introduction. Our uniqueness proofs are based on the following 
orthogonality relation involving solutions in Q and impedance maps on dQ. 

Proposition 3.1. Let (/ii,7i) and (/i2,72) be two pairs of coefficients such thatuo is not a 
resonant frequency for either. Let Ai and A2 be the corresponding impedance maps. Then, 
for any Xj =( if* |0 i?j ) satisfying 

iP^iD) + V,„,,)X, = Q 
(P^(D) + V;„-JX2 = 

in Q one has 

[ X;{V,,,,, - V,,,,,)X, dx = i f (A2 - Ai)(z/ A ifi) ■ H2 dS. 
Jn Jan 
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Proof: If 

{u\v) = / v*udx and {u\v)Qn = / v*udS 
Jn Jon 

the integration by parts identity 

{P^{D)X\X') = -i{P^{iy)X\X')9n + {X\P^{D)X') 

holds. 

Since *j = = for j = 1, 2 one has that {Xi\V*^ ,^^X2) = {Xi\Vn2^^^X2), hence 

- v,,,^,)x,\X2) = (v;,,,Xi|X2) - (Xi|\/;,^^X2) 

= -(P^(D)Xi|X2) + {X,\P^{D)X2) 
= t{P^{iy)Xi\X2)dn 

= i [ E2- {ly AHi) -H2- {ly A Ei) dS 

JdQ 

Let X*2= [o HI\q ^2) 

be the solution of the boundary value problem 

{P^{D) + V^„^,)X2^0mn 
1/ A H2 — ly A Hi on dft. 

Since X2 on the boundary verifies 

E2-{uA Hi) ^E2-{uA H2) 

we have 

{{V^,a. - ^^2,72)^11^2) =t ! E2-{yAH2)-H2-{yA ^1) dS. 

JdQ. 

Finally, note that 

-i{P^{u)X2\X2)xi = {P^{D)X2\X2) - {X2\P^{D)X2) = -{V^,,^,X2\X2) + {X2K,,^,X2) 
= -(^M2,72^2|-^2) + (-'^2|V'^*2,-y2-'^2) = ((^M2,72 " ^M2 ,72 ) ^2 1 -^^2 ) = 0. 

Here we used again that ^'2 = $2 = ^2 = ^2 = 0. Prom the last identity one gets 

i I E2-{vAH2)dS^i I H2-{iyA E2) dS. 

JdQ. JdQ. 



Therefore 



((^Ml,7l - ■^M2,72)^l|^2) I H2-{VAE2)-H2-{VA E,) dS 

Jan 

= i [ H2-{A2-Ai){u A Hi) dS, 
Jan 



Ian 
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which is the desired integral formula. 



□ 



We next plug the solutions constructed in Proposition 12.31 into the integral formula given 
in Proposition 13.11 Since Aia = on F for a E TH^^'^{dVt) with supp(a) C F, we have 

/ ^l{y^^ull - ^M2,72)^i dx = i f (A2 - A ■ ^2 dS 

Jn Jdn 

= 1 [ (A2-Ai)(i^Af)i)-^2C^5. 

J To 

Note that A i32|ro = means that the tangential components of Sj2 on Fq vanish. On the 
other hand, (A2 — Ai)(z/ A Sji) only has tangential components. Hence 

(A2-Ai)(z/Aioi)-:^2 = 

on Fq. Since the boundary term in the integral formula vanishes we get the relation 

/ •^2(^1,71 ^(12,12)^^ dx = 0. 
Jn 

Following the same notation we write 

Yj{x) = dia.g{— I 4, 14) Yj{x{x)) 

and we thus have 

Xi = Xi + Xi = diag(/i-'/', 7r'^')(>l + Yi) 
X2 = X2 + X2 = diag(/i2"'/',72"'^')(>'2 + Y2). 

Since for r large enough the scalar components vanish, i.e. X* = Xj = for j = 1, 2, 
writing Vf,,^^, - V^^ as in (11.131) gives 

^2(^m,7i - '^/^2,72)^i = X;diag(cc;(/ii - fi2),^{li - 72))Xi 

= (F2 + 5>2)*diag(/I,7)(n + yi) 

where 

/il - /i2 ~ 71-72 

If now we go back to the integral relation, we get 

0= / Xl{V^,,,,-V^,,,,)X^dx= [{Y2 + Y2ydmg{fl,^){Y,+Y,)dx 
Jn Jn 

= / Y2dia.g{Jl,'j)Yidx + / Y2^iag{'jl,'j)Yidx + / Y2dia:g{Jl,'y)Yidx 
Jn Jn Jn 

+ / F2*diag(/x, 7)Y'i (ix. 
Jn 

18 



We gather the first two terms and the last two terms (cross terms) in different expressions. 
Observe that 

y'2*diag(/1, 7)11 dx= {Y2{x))* diag{— 14, J4)diag(/I, 7)diag(— /4, 14)^1 (x) dx 
Jn 

(F2(a;))*diag(/I, ^)Yi{x)dx 
F2*diag(/1, 7) Yi dx 

hence 



n 



y2*diag(/i, 7)Yi (ix + / Y2diag{^,'y)Yidx = / Y2diaLg{^,'j)Yi dx. 
Jn Jo 

With a similar argument the cross terms become 

F2*diag(/I, 7)Yi + / y2*diag(/x, 7)Y'i (ix = / F2*diag(/x, 7)Yi (ix. 
Jn Jo 



Therefore, 



0= / X;{V^,,^,-V^,,^,)X,dx= [ Y*dmg{Jl,Widx+ [ Y*dmg{]l,Widx. (3.27) 
Jn Jo Jo 

Below we will integrate by parts in O many times. Even though O may not be a 
Lipshitz domain, this procedure is easily justified by integrating by parts separately in the 
C^'^ domains fl and Cl and by noting that the contribution from {X3 = 0} vanishes. Now, 
the terms on the right hand side of (13.271) may be simplified by writing 

/ {Y* + Y*)diag{]l,Widx= [ {Y* + Y,*)dmg{Jl,^){P^{D) -W'^^^^jZ.dx. (3.28) 
Jo Jo 

Note now that 

P^(D)(diag(7,/I)^i) = dmg{Jl,^)P^{D)Z, + {P^{D)dmg{^,Jl))Zi 

and 

/ (r; + r;)^^(^)(diag(7, Ji)z^) dx = 

Jo 

= [ {P^{D){Y2 + r2))*diag(7,/I)Zi dx-i ! {Y; + r2*)^^('^)diag(7, aI)^i dS 

Jo JdO 

integrating by parts. We know by (10.41) that yUi = ^2 and 71 = 72 on F, and this is valid 
on do as well. Since Y2 + Y2 is solution for the rescaled system 



[P^[D) + W,,^^MY2 + Y2)=Q 
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in O, we have 

/ (r; + F;)P^p)(diag(7,/I)Zi)rfx = - / iW,,,^^{Y2 + Y2)rdmg{^,Jl)Z,dx. 
Jo Jo 

From f l3.28p and the above observations we get 

/ {Y* + r;)diag(/I, Wi dx = (3.29) 
Jo 

= - [ (r; + F;)(iy;^,^^diag(7,/I) + (P^(D)diag(7,/I)) + diag(/I,7)W^;,.,J^icia;. 
Jo 

It is a straightforward computation to check 

W;^^^dmg{^,]l) + (P^(D)diag(7,/I)) + diag(/I,7)lV;„,, = 
K2diag(7, Jl) + /«idiag(/1, 7) + P^(Z^7, ^/^) + ^^(^7, DJl), 

where 

A^i + ^ 71+72 



(/^l/i2)l/2' ' (7172)1/2 ■ 

From now on we denote for short 

U := K2diag(7,/I) + Kidiag(/I,7) + P^(P'7, D^i) + P^{D^, DJi). 
Observe that these computations could be written separately for Y2 and 1^2, giving 



y2*diag(/I,7)yic/a; = - / Y^UZ^dx, / r2*diag(/1, 7)^1 da; = - / Y^UZ^dx. 
o Jo Jo Jo 

Recall that solutions of the rescaled Maxwell system and the Schrodinger system look 
respectively like 

F2 = e'^^'-^Yl + Y^ + y;2), Zi = e'«^-"(Zi + Z\ + Z]). 
Using the estimates in Proposition 12.11 and in 12.181 we have 

lim / K,*diag(/I,7)yirfx = - lim / e'^^^-~^^^ ''{Y^YUZldx 

- lim / e'^^'^'^^^-^Y^yUZldx 

- lim / e'^^^-~^^^%Y^YUZ\dx. 
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Noting that Ci ~ C2 = ^ taking the hmit in fl3.27p . we obtain 

= - hm / e'^ ^'iY^yUZ^dx 
^^+°° Jo 

- hm / e'^ %Y^yUZ^dx 

- hm / e'^ '^iV^yUZ'^^dx 



T— > + 00 



- hm / Y;UZidx. 



(3.30) 
(3.31) 
(3.32) 
(3.33) 



o 



4. Some technical computations. In this technical section we compute the limit of 
the terms given in dS^OD, dSSID, and flX^ . 



Lemma 4.1. The limit in ^3. 30\) is given by 

lim / e'^ ''(YlyUZ^dx = 
= fc((C-62)(C-&i) + (C-«2)(C-ai)) / e'^-^{K^ + K2)(ii + ^)dx 



o 



+u;{(.hyC-h) I e^«-(-A)(-)'^'^a; + c.(C-a2)(C-ai) [ e^^-^{-A)(^Y' dx 

V/i2/ Jo V72/ 



O 



+2kuj j e 
'o 

+2ku I e'^-^ 
'o 



-{(■ h){CA ai) -Df^y^' + (C ■ a2)(C A h) ■ d(^) 



-(C • &i)(C Aa,).D (^Y' + (C ■ aO(C A h) ■ d(^) 

\fi2/ V72/ 



1/2 



1/2] 

rfa;. (4.34) 



Proof: Remember that 

with 
and 



m 



* = (0 6*|0 a\) 



A simple commutation in the middle term after writing the definion of U , gives 



ei^--{Y^yUZldx=- I e'^-%P^iC2)Z^yUmidx 
T Jo 



(4.35) 



+ 



2r 



e^i-i^ZlyP^iC^yP^iCi) + P^(Ci))(/t2diag(Al,7) + «:idiag(7, y5))mi dx (4.36) 



2r 



e^^-\ZlyP^{C2){P^{D^,Dri) + P^{D^,DJi)){P^{Ci) + P^{Ci))m,dx. (4.37) 



3± 
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Before computing the limit of fl4.35p let us write 



( k(2 ■ Ci2 \ 

(C2 • h)C,2 - kC2 A a2 
K2 ■ &2 

V(C2-a2)C2 + K2 A62/ 



(4.38) 



It follows from this expression that 

lim - / e'^-''{P^{C2)ZlYUm^dx = 

r~.+oo T Jo 

= fc^e^^-"(o (C-&2)Cl0 (C-a2)C*) (/t2diag(7,/I) + /€idiag(/I,7)) 
+k e^«-^ ( (C ■ 62)Cl0 (C ■ a2)C ) (P^iDl, m + P^{D^, DJl)) 



/0\ 

h 


h 




dx 



dx 



k I e'«-" 
'o 



(C ■ ^2)(C ■ &i)(/«27 + + (C ■ a2)(C ■ ai)(/«2/i + fi:i7) 

1/2 

./ii/ - \7i 

To compute the limit of the term in (14. 361) observe that one has the relation 



Jo L ^^1^ ^71 



dx 

dx. (4.39) 



^(^o')*^^(C2)(^^(Cl) +P^(Cl)) = ^2^*2P^C2)iPHCl)+P^iCl)) 

+ ^m;{P^(C2) + P^(C2))P^(C2){PHCi) + P^iCi)). 



(4.40) 



Using the commutation formulas fll.91) . fll.lOl) and (11.111) . and the fact that Ci ~ C2 = 'C? 
can see that 

_ (^^^((2) + ^^(C2))^_^(C2)(^^(Cl) + P^(Ci)) = 

= (Ci ■ Ci)P^(C2) - P^{0P^{Ci)P^{C2) + (Ci ■ Ci)i'^(C2) - P^iOP^iQP^iCi) 

+ (C2-C2)(^^(Cl)+P^(Cl)) 

= -P^iOP^iOP^iQ - P^iOPHQP^iCi) + 0{r). (4.41) 
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Having in mind formulas (I4.4UI) and fl4.4ip . we compute the limit of (14.361) : 

lim 1 / e^«-(Z2)*P^(C2)(P^(Ci) + i'^(Ci))(/^2diag(/I,7) + «:idiag(7,/I))micix 

r^+oo It Jq 

= ^ I e*«-^m;P^(C)(P^(C) + i'^(C))('«2diag(/I,7) + /«idiag(7,/I))midx 



Jo 

e*«-^m;(P±(0 + P^(0)P^(C)i^^(C)('t2diag(/I,7) + «:idiag(7, /I))mi dx 



o 



_ k 
~ 2 

k I e'^-^ 
o 



e*^''m;P^(C)P'^(C)('^2diag(/i,7) + Kidiag(7, /i))mi dx 



o 



(C ■ &2)(C ■ &i)(/«2Ai + i^il) + (C ■ a2)(C • ai)(/«27 + /^i/") 



(4.42) 



Here we used P^(C)P^(C) = C ■ C = 0, 



p±(C)p^(c) = p^iOPHc) = 2 



/ 



V 













oc®C/ 



and 



(p±(0 + P^(0)^^(C)P^(C) = o. 

Finally we study the limit of (14.371) . Let us write 



where 



We have 



P^(e2(r))^-P^(0asr^+oo. 



2r 



- / e'^-%Z^yP^iC,){PHD^,Dfl) + P^{D^,DmPHCi) + P^{Ci))midx 



o 



At 



i- I e^^%ZlyP^{Ur)){P^{Di,Djl) + P^{D^,DmP^{Ci)+P^{Ci))^idx 



o 



+^ / e'^-^{ZlrP^{0{P^{D^,Dri) + P^{D^,DJi)){P^{Ci) + P^{Ci))m,dx. 



o 



We split Zg into pieces depending of order of r, 



^0 



/ C2 ■ ^2 \ 



C2 ■ &2 

V / 



1 

+ - 

r 



/ \ 



\ka2 J 



:= ^M^ + -Ml 

T T 



(4.43) 
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and we set 



/C-«2\ 



Ms = lim 

T—*+00 







V J 



It follows that 
lim — 



e^^-^iZ'^yP^^iQiP^Dj, DJI) + P^{D^, L'/I))(P±(Ci) + P=^(Ci))mi dx = 



o 



e'^-^M*P^{0{PHDl, D^) + P^(D7, L'A^))(P±(C) + ^^(0)^1 



Observe that for any a e the formula 

P^{a){P^{D^,Dti) + P^{D^,DJi)){PHC} + P^{0)mr = 
/ 0-^(71/72)^/2 
-a A D(-fi/^2y/'' 



4a; 



v 















aAD(/Xi//X2)V2 oy 



V y 



holds. Using that C2 = '^C ~ |C + C'(r the limit is 



lim — 

T— >+00 2t 



o 



+2ku) I e 
'o 



'7iy/2' 
^72/ 



dx 
dx 



-u ! e^^- (C-&2)(C-M(e-^)(-)'^' + (C-a2)(C-aO(e-i^)(^) 
Jo L ^/^2/ ^72^ 

-(C • 6i)(C A 02) • D (^) + (C • aO(C A 62) • ( 
+ hm — / e^«-4(C2-62)(Ci-&i)(^i^)(-)'^V(C2-a2)(Ci-ai)(^i^)(^)'^' 

T Jo L ^A*2^ ^72^ 

In the first term on the right, integration by parts gives 



dx 



-u; / e^«- [(C-62)(C-M(e-^)(-)''" + (C-a2)(C-ai)(e-D)(^)'^' 
Jo L ^/^2>' ^72^ 

= / e^«- [(( • 62)(C • 6i)(-A) (^) '^V (C • a2)(C • ai)(-A) (^) 
Jo L ^/^2/ ^72^ 

(C ■ &2)(C ■ • D) (^)"' + (C ■ a2)(C ■ ai)(z. ■ (^)'^' 

\U,r>J ^72^ 



+za; / e 

'dO 



dx 
dx 
dS. 
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Since dlfii = dlfi2 and (9^7i = 9^72 for / = 0, 1 on F, the boundary term vanishes. On the 
other hand, for the last term on the right we integrate by parts to obtain 



o 



dx 



o 







Jo 


_V/i2/ 



dx 
dx = 



and 



72^ 



1/2 



o 



72^ 



1/2 



dx = 
dx = 0. 



We introduced —1 in order to get rid of the boundary terms in the integrations by parts. 
Eventually one has 



o 



uj e 
'o 



(C-&2)(C-&i)(-A) 



/iiy/2 



(C-a2)(C-ai)(-A) 



7ixi/2' 



72/ 



+2ku [ e^«- L(C.6,)(CAa2)-/^(^)'^' + (C-ai)(CA&2)-/^(-)'^' 

Jo L ^/^2^ V72/ 

Putting (14.391) . fl4.42p and fl4.44p together we get the final result. 
Lemma 4.2. The limit in ^3. 31\) is given by 

lim / e'^ %Y^yUZldx = 
(C ■ a2)(C ■ ai)(/«27 + + (C ■ 52)(C ■ &i)(t2Ai + i^il) 

72 V72/ /i2 V/i2^ 



rfx. (4.44) 
□ 



'o 



(C ■ a2)(C • fli) 



1/2 



da; 



+ e 
'o 



(C ■ if )(C ■ ai)(«:27 + /^ipf) + (C ■ i?2'')(C ■ hl){K2{l + /Cl7) 



dx 



+2u [ e^«- [(C ■ b,){R^ A • /^(^)'^' - (C ■ A f) ■ 

Jo L ^/^2^ V72/ 

-o;(C-6i) [ e'«-"2(C ■ D)i?| f ^V^' - 1 dx 



dx 



-a;(C ■ ai) / e*«-^2(C ■ /^)i?| 



dx, 



(4.45) 



25 



where 



with 



K2 



R2 = lim tZ_^. 

r— >+oo 



Proof: Since = P^{C2)Z^^ - Wj^^^^^Z^ we have 



o 



e^i-^iY^yUZldx = / e'^-''{Zl^YP^{C2)UZldx 



o 



- [ e'^-^{zlyw,,,^^uzldx. 

Jo 



We can write 



and set 



Ml = lim Zl 





V / 



M2 = lim 



/C-a2\ 



C-b2 

V / 



Then the hmit of flOTD is 



lim - [ e'^-%Zl)*W^,,^^UZldx = 

= - [ e'^ ''K2M;UMidx+l [ e'^ ''M;P^{DP2,Da2)UMidx 
Jo 2 Jq 

= - e*^'^K2M2('^2diag(7,/i) + Kidiag{Jl,'^))Midx 
Jo 

1 [ e'^-''M;P^{Dp2,Da2)iP^iD^,Dil) + P^{D^,DJl))Midx 

2 Jo 

(C ■ a2)(C ■ ai){K-2l + i^ipj + (C ■ 52)(C ■ bi){i^2l^ + ft^iT) 



(4.46) 
(4.47) 



o 



+00 j e 
'o 



(C ■ ^.)(C ■ ■ + (C ■ m- h)^ ■ d(^) 

72 V72/ yU2 V;U2/ 



1/2- 



dx 

dx. (4.48) 



In the second equality we just wrote down U and we used 

M;{P^{D^, Dfl) + P^iD^, DJi))Mi = 

and 

M;P^{DP2, Da2)iK2diag{^,ll) + Kidiag(/1, 7))Mi = 0. 
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We now compute the limit of the right hand side of (I4.46P : 

hm / e'^-^{Z\yP^{l^)UZldx = 

Jo ^ 

+2UJ [ e*«- [(C ■ A -dI^Y' - (C ■ ai)(i?? A() ■ dI^)' 

Jo L ^/"2^ 

+2UJ I e^«- [(C-&i)^|(C-/^)(^)'^' + (C-ai)^|(C-/^)(^)'^' 
Jo V ^At2^ 

Integrating by parts one gets 

/ e^^-E:i{C-D){^)"" dx= [ e^^-R^C-D) \(^) 

Jo ^At2/ Jo L^/"2/ 

= - / e^«-^(C ■ D)R^ 
Jo 



dx 



dx 



dx. 



dx 



and 



! e'^-^R!(X-D)(^Y'dx= [ e^^-^RtiC-D) 
Jo ^72^ Jo 

e'^%C-D)R 



o 



72^ 



dx 



dx 



dx. 



Once again, we introduced —1 in order to get rid of the boundary terms. Finally 

lim / e'^-''{Z\YP^(X^)UZldx = 

= f e^«-" [(^^f)(^al)(/^27 + '^l/^) + (^^2'')(^&l)('^2/^ + /^l7^ dx 
Jo L 

+2C0 [ e^«- \(C-h)iRiAC)-D(^Y'-{C-ai)iR^AC)-D(^y 

Jo I ^1^2^ 

-2uj(C-bi) f e*«"(C ■ f^V^'-l dx 

Jo L^/^2^ 

-2uj(C-ai) f e'^^'iC- D)Rt f^V^'-l dx. 

Jo L^72^ 

Putting (14.481) and fl4.49l) together we get the final formula. 



dx 



(4.49) 
□ 
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Lemma 4.3. The limit in ^3. 32\) is given by 

lim / e'^-''{Y^yUZ\dx = 



O 



dx 



+2UJ I e^«- [(C.a2)(CA^f)-Dp)'^'-(C-62)(CAi?f)-/^(^) 
Jo L ^7i^ ^/^i^ 



1/2- 



dx 



'o 





J 1/2 








1/2 


.V72^ 





dx 
dx, 



(4.50) 



where 
with 



R\ = ( i?f (i?f)* iRf (i?f^* 



i?i = lim rZ\^. 

r— >+oo 



Proof: We use the computations done in fl4.38p for = P^{(2)Zq and look at the limit 
lim / e'^-''(Y^yUZ^_^dx = lim / e'^-%P^{(2)Z^yUZ\dx 

= J^^'^'" (o {C-b2)C\0 (C-a2)C) URidx. (4.51) 
We insert the definition of U in (14.511) and we compute 



O 



(C-&2)Cl0 iC-a2)C)URidx 



o 



+ /^e*^-^ (^0 (C-&2)ClO (C-a2)C*j (P^(/^7,^/^) + i^^(^7,^/^))^irfa; 

dx. 



{^21 + /tl/i)(C ■ &2)(C ■ Rl) + ('«2/i + /€l7)(C ■ a2)(C ■ 



dx 



[(«:27 + «^i/^) (C ■ 62) (C • O + ('^2/2 + «:i7) (C ■ ^2) (C ■ ) 

/ e'«- fef(C-&2)(C-/^)(^)'^'-(C-&2)(CA^f)-I^P)'^' 

/^aN 1/2- 



+2cj / e^«-^ 
'o 



RfiC- a2)(C- D) (^) + (C ■ a2)(C A i?f ) ■ . . 
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Integrating by parts one gets 

Jo ^1^2^ Jo 



and 



Jo ^72^ Jo 1^12^ 



- / e^«-(C ■ D)Rt 
'o 



dx. 



This gives the final result. 

The last technical lemma shows that the cross terms do not contribute to the limit. 



□ 



Lemma 4.4. The limit in ^3. 33\) satisfies 



lim / Y*UZidx = 0. 



(4.52) 



Proof: Recall that 

Y* = e-<2-(i;2 ^ y2 ^ Y2y 2^ ^ e^'-\Zl + Z\ + Z]), 

where Yj = diag(— /4, l4^)Y^{x) for j = 1, 0, r. Note that writing the dot product in terms 
of the basis / in f l2.26p one has 

Since we chose /s = 63 we have that x^ = X3, hence 

/ Y*UZidx= [ e"t'-{Y^ + Y^ + Y^yU{Zl + Z\ + Zl)dx 
Jo Jo 

2(^2 ^ P)i/2Mj"' ^ D^^e'^^iY^ + Y^ + Y^rU{Zl + Z\ + Z]) dx 
2(r2 + ' e^'^^D^.liY^ + Y^ + Y^YU^Zl + Z\ + Z].)] dx 

^ \^\ ^ Jao 
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The boundary term vanishes because we assumed that d^Hi = dlfi2 and (9^7i = 9^72 for 
/ = 0, 1 on r. Since and Zq are constant vectors with respect to x we have 

D,,[{Y,' + + VrU{Zl + Z\ + Zl)] = D^,{Y^ + VTU^ZI + Z\ + Zl) 

+ (Y^ + Y^ + Y^rD,,U{Zl + + Zl) 

+ {Y^ + Fo' + Y^yUD,,{Z\ + Z^). 

Note that from the asymptotic behavior in fl2.18p and in Proposition 12.11 one can see that 

d,,{yS + Y^yuizl + z\ + zl) = D.,,{Y^yuzl + 
{Y^ + + i;')*^x3^(^o' + z\ + z^) = {YlyD,,^uzi + 

(1>2 ^ ^ Y^)*UD,.^{Z\ + = {Y^yUD^.X + 
in L'^{0); moreover, 



D,,iY,')*UZl, 



LHO) 



0[t 



{jiyuD,,zi 



L\0) 



0[r 



On the other hand, we know by f l4.38p and f l4.43p that 



2\t 



( (C2 • &2)a|0 (C2 • a2)d ) + 0{1) := -{LlY + 



in LF'{0). Therefore 



r 



( Ci ■ ai 0|Ci ■ 61 ) + 0{t-') := -(M^i)* + 0{r-^) 



hm / 



o 



hm (2(r^ + A;2)V2|^)"' 
-^+00 V 1^1 / 



e"^^D^^{Y^yUZldx 

+ /" e'^^{Y^yD^^UZldx+ [ (^^^{Y^yUD^^Zldx 
Jo Jo 



hm 2(r" 

T— > + 00 



- hm -[2{T' + k'y/' 



1^1 



e^<^-(ri'(i;))*diag(-/4, ii)D^.^UZl dx 
e*'^^(L^)Miag(-/4, iA)D^.^UMl dx. 



o 



Note that 



{Llydi^g{-h,h)D^,UMl 



2uj{C, ■ a2)(Ci ■ «i)(C2 • D)D,; i-V' - MC2 ■ &2)(Ci ■ h){C, ■ D)D,, 



y/2 



Since the Riemann-Lebesgue lemma imphes that 



hm 

r— ++00 



o 



/i2 



hm 

T— > + 00 



o 



7i\ 
72^ 



1/2 



(ix = 0, 



we have proved the resuh. 



□ 



30 



5. Proof of Theorem 1. Assuming the conditions in Theorem 1, we have proved that 
the identity in fl3.30p - fl3.33p holds and that the terms in that identity have the hmits given 
in Lemmas 14 . 1 1 to In this section we show the equations verified by the coefficients and 
then give a proof of Theorem 1. 

Proposition 5.1. Let ^ g M^. With appropriate choices of the constant vectors ai, 02, 
bi and 62 one has the following equations: 

(a) If we choose 61 = 62 = C o,i^d, ai = 112 = ( 



o 



dx = 0. 



(5.53) 



(b) If we choose ai = 02 = C ^'^^ ^1 

'1 



o 



dx = 0. 



(5.54) 



Proof: In order to get these equations from the previous lemmas we will need the equations 
f l2.24p and fl2.25p . Let us write these equations by components: 



2{C-D)Rf = 

2{C-D)Rf = 
[2{C-D)hRf 

2(C ■ D)Rl = 
2(C ■ D)hR^ 
2(C ■ D)R^ = 
[ 2(C ■ D)hR^ 



|Aai - Kl + k'^ - liDai] 
2(C ■bi)DK^ 

(|A/5i-/€? + A;2-i(D/?i) 
2(C-ai)D/€i, 



iAa2 -nl + k"^ 



2(C ■ b2)DK2 

■(iA/52-7€2 + P 
2(C ■ a2)DK2. 



\{Da2? 



(C ■ ai) 

(C • a2) 
(C ■ b2) 



(5.55) 



(5.56) 



We take complex conjugates in fl5.56p in order to get the information required. 

2(C ■ D)Rl = (lAa2 -Kl + k^- \{Da2f) (C ■ ^2) 



2(C ■ D)hR^ = 2(C ■ b2)DK2 

2(C ■ D)R:I = f iA/32 ~^l + k'- UDf32r) (C ■ b: 



(5.57) 



[2{C ■ D)hRi = 2{C-a2)DK2. 



These equations are valid in 
and fl5.57p by one gets 



^ Multiplying the second and fourth equations of f l5.55p 



C-i?f = (C-6i)(/s:i-fc) 
C.Rf = (C-ai){K,~k) 
C-R^ = {C-b2){K2-k) 
[C-R^ = {l-a2){K2-k). 



(5.5J 
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To see the last identities we used the decay of Rj for j = 1,2 at infinity, and the fact that 
( ■ Rf — (C ■ C ■ Rf ~ (C ■ C ■ R2 ~ (C ■ b2)K,2 and C ■ i?f — ■ 112)^2 are entire 

functions in the variable C, ' ^■ 
Let us denote 



(Tl) = 




hi) 


lo 


-2(C- 




(T2) = 




ai) 


[ e' 
Jo 




D)Rt 


(T3) = 




h) 


[ e'' 
10 


■-2(C 


D)Rf 


(T4) = 




02) 


1 e' 
Jo 


^^■2(C 


D)Rf 



fJ'2^ 

71 X 1/2 

fJ'2^ 
71 X 1/2 

72^ 



fix 
da; 



which are terms from (14.451) and f l4.50p . Then from the first and third equations in f l5.55p 
and f l5.57p one gets 



(Tl) + (T2) + (T3) + (T4) = 



o 

-t^(C ■ ai)(C • 02) / e 

JO 



V/i2/ 



Noting that 



2 Vo'i/ Vo''j/ Vo'i/ V72/ 



one can write 



(Tl) + (T2) + (T3) + (r4) = (TS) + (T6) 



-.(c.M(c-.2)Xe-[(;^)%(;^)-..( 



1/2 



:A(/32-/5i) 



where 



(T5) = -^(C-fei)(C-&2) / e 



'^2 + '^?-^(W + ^PAf 



(T6) = -uj(C- ai)(C ■ 02) / e'«-" - + _ l(/)a2)' + ^pai) 

'o 



4 



^ ~ ^A(a2 -ai) 
V72/ 2 























.^72^ 





dx 
dx, 

dx 
dx. 
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Note that 



72^ 



1/2 



1^2 

dCA''' .D(^-lf" = \Da2 ■ Df^V'' - \Da, ■ D^A 



y/2 



7^^/32 



1/2 



V/i2/ 

(-) 

V72/ 



^2^ 



1/2 



V72/ 



Denoting 



(T7)=^(C-6i)(C-&2) / e^«-"^A(/52-/3i)cia; 



we get 



(T8) =cu(C-ai)(C-a2) / e'«-^'^A(a2 - «i) c^a; 

Jo ^ 



(Tl) + (T2) + (T3) + (T4) = (T5) + (T6) + (T7) + (T8) 



V/i2/ 



V2/ 



-^(C-&i)(C-52) / e^«-^' 
-.(C ■ aO (C ■ a,) e'^- [ifla, . (^) "''^ - . D (2l) 



1/2 



da; 



-c.(C-&i)(C-fe2) / e*«-^(-A)(^)'^'rfx 



io ^72^ 



We consider the next terms from fl4.34p and (14.451) 



(T9)=a;(C-62)(C-&i) / e'^\-^)(^^^'^ dx 



o 



(T10) = a;(C-a2)(C-ai) 



Jo ^72^ 



(Til) = uj I e*«-" 



(C-a2)(C-ai)^«2-/^ 



/7i\i/2 
V72^ 



fix. 



hence 



(Tl) + (T2) + (T3) + (T4) + (T9) + (TIO) + (Til) = (T5) + (T6) + (T7) + (T8) 

+ iBft.fl(^ii)"^^ 

2 V/i2/ 



+C.(C-6l)(C-&2) 

+a;(C ■ ai)(C ■ ^2) 
We note that 









^/^\ 


Jo 






^/i2/ 


f gif- 


-Da;2 


^( 


^7i\ 


Jo 


.2 




^72^ 



1/2 



dx 
dx. 



1/2 



1/2 



2 V72y 



{Dai — Da2) 
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hence 



(Tl) + {T2) + (TS) + (T4) + (T9) + (TIO) + (Til) = (T5) + (T6) + (TT) + (TS) 

Jo 4V/i2/ 



o 



(5.59) 



Considering the terms 



(T12) = / e 



(C ■ )(C ■ ai)(/t27 + /^i/^) + (C ■ R2)iC ■ bi){K2jl + na) 
(«:27 + «i/^)(C ■ &2)(C • Ri) + {f^2ji + Ki7)(C ■ a2)(C ■ i?f ) 
(T14) = A;((C-62)(C-6i) + (C-a2)(C-ai)) / e*«-^(/s:i + /t2)(/I + 7) 



(T13) = / e 



dx 
dx 



o 



(ri5) = - / e'«"fi:2 (C ■ a2)(C ■ ai)(fi;27 + /^i/u) + (C " &2)(C • bi){n2f^ + ^17) c/x 
Jo -I 

from (14.341) . (I4.45P and (14.501) : and using (15.581) we have that 

(T12) + (T13) = 

(C ■ «2)(C ■ ai)(«2 - k){K2^ + KiJl) + (C ■^2)(C ■ ^l)(«2 " k){K2jl + K^^) 

(K27 + Ki/i)(C • 62)(C ■ &i)('ti - A;) + {K2JI+ K,rf)(C-a2)(C- ai){Ki - /c) 



J^-x 



o 



'o 



dx 
dx. 



Therefore 

(T12) + (T13) + (T14) + (T15) = 
= (C ■ a2)(C ■ oi) / e*^'^Ki(fi;2Ai + K17) c/x + (C ■ &2)(C • &i) / e^^'"" ki{k2^ + kiJi) dx. (5.60) 



Now we make the choice 61 = 62 = C ^"^^ ai = 0,2 = (■ The equations (I5.55P and (I5.57j) 
imply that Rj" = for j = 1,2. Considering all the terms from the technical lemmas in 
Section HI we see that the only surviving terms in (I3.30I) - (I3.33I) are those which appear on 
the left hand sides of (15.591) and (15.601) . Thus these computations give 

/ e^«- \\a{P2-Pi) + ](-y^\iDp,f - {DP2Y] + ^{n2^ + n,Jl) 

Jo ^ ^^^2^ ^ 

dx = 0. 



1^2 
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If we choose Oi = 02 = C 61 = 62 = C) "we get that = for j = 1, 2 for the same 
reason as above. We consider again all the terms from the lemmas in Section H] and we get 

4 V72- 



O 



^A(a2 - «i) + 7 (z^y^^iDa.y - {Da^Y] + ^{k^JI + k^) 



+ (4 



1 



-pas) - -{Da 



4^ 4- 
Finally, by a direct computation we see that 



72' 



1 



dx = 0. 



1/2 



1/2 



0, 
0, 



hence the result is obtained. 

The preceding arguments are valid for all ^ G with |,^'| > 0, since this was the 
assumption on ^ in the construction of reflected CGO solutions. However, the continuity 
of the Fourier transform proves the result for all ^. □ 

The above proposition shows that the coefficients satisfy the following equations in O: 

-|A(ai - 02) + jD{ai + 0:2) ■ D{ai - ^2) + t^^(/Wi7i - /W272) = 0, 
-iA(A - P2) + lD{f3, + P2) ■ D{/3i - P^) + ^'(/ii7i - m2) = 0. 

Let u = (71/72)"'^^^ and v = (/ii//i2)^^^- The equations become 

-A(log u) + (7i72)"^''^^(7i72)^''^ ■ ^(log u) + uj'^ifii'ji - ^1212) = 0, 

-A(log V) + (/ii/i2)-'/'/^(/il/X2)l/' ■ D{log V) + UJ^{lli-ii - /i272) = 0. 

Multiplying the equations by (7172)^^^ and {fj.iH2y^'^ we obtain 



D ■ ((7i72)'/'^(log u)) + uj\fx^^^ - /i272)(7i72)'/' ■ 

D ■ {{fiil22y/^D{\og V)) + Uj\fii-fi - 1^272) (;Ui/i2)^/^ 



0, 

: 0. 



Finally, these may be written as 



D ■ {-12DU) + cu2/i27i(w'^' - l)w = 0, 
D ■ ifX2Dv) + uj'^iil-i2{u^v'^ -l)v = Q 

in O. This is a semilinear elliptic system for u,v E C^(0). The assumptions in Theorem 
1 imply that u\do = v\do = 1 and duu\so = duvlgo = 0. The following formulation of the 
unique continuation property will therefore imply that 71 = 72 and /ii = fi2 in O, thus 
proving Theorem 1. 
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Lemma 5.2. Let a, b be non-vanishing complex functions in C^{0) with positive real part, 
and let p,q be complex functions in L°°{0). Suppose that u,v E C^(0) satisfy in O 

D ■ {aDu) + p{u^v'^ - l)u = 0, 
D ■ (bDv) + q{u'^v^ - l)v = 0, 

and also u\do = ^|ao = 1 o,nd duu\do = duv\so = 0. Then u = v = 1 in O. 

Proof: Let u = 1 + a^^l'^z and f = 1 + b^^l'^w. Using the identities 

D . {aD{a-y^z)) = a'/\-Az + ^z), 

D ■ ibDib-^'^w)) = b'/\-Aw + ^^w), 
u'^v'^ - 1 = {aby^/'^{zw + b^/^z + a^/^w){uv + 1), 

we obtain the equations 

-Az + ^^z + p{zw + b"^z + a^/^w) = 0, 

-Aw + + q{zw + b^l^z + a^/^w) = 0. 

Here p,q & L°°{0) are given by 

p = a~^/'^p{ab)~^/'^ (uv + 

q = b~^/'^q{ab)~^/'^{uv + l)v. 

Thus, there is C > such that the following differential inequalities for z^w & C^(0) are 
valid almost everywhere in O: 

\Az{x)\ < C(\z{x)\ + \w{x)\), 
\Aw{x)\ < C{\z{x)\ + \w{x)\). 

We also have -zjao = "U^lao = and dyz\Qo = d„w\Qo = 0. The unique continuation property 
holds in this setting, which can be seen by applying a scalar Carleman estimate to both z 
and w (for details see |KSaUj ). Thus z and w vanish identically in O. □ 

6. Proof of Theorem 2. Here we prove an analog of Theorem 1 for the case where 
the inaccessible part is part of a sphere. Following the idea in [I], this reduces to the 
hyperplane case by using the Kelvin transform. 

We first need to check how the Maxwell equations behave under the Kelvin transform. 
For this it is convenient to consider E and H as complex valued 1-forms in (that is, a 
vector field (i?i, £^2, E^) is identified with the 1-form Ei dx^ + E2 dx"^ + E3 dx^). We equip 
Q with the Euclidean metric e, and introduce the Hodge star operator *e which maps a 
1-form Xi dx^ + X2 dx^ + X3 dx^ to the 2-form Xi dx"^ A dx^ + X2 dx^ A dx^ + X3 dx'^ A dx'^. 
If d is the exterior derivative, the Maxwell equations in form notation are 

{dE = iujii *e H . „ 
dH = —iLwy *e E 
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Assume that is not in Q. The Kelvin transform is given by 

x{x) — \x\~'^x, F{x) — x{x) — \x\~'^x. 
If Q = F'^{Q), then F is a conformal transformation from {Q, e) onto (^2, e): 



F*e = \x\ 



-4, 



Let E = F*E, H = F*H, ji = and 7 = F*7. The following is the transformation law 
for the Maxwell equations under the Kelvin transform. 

Lemma 6.1. One has 

dE — iuja *e H . „ I dE = iujjl\x\~'^ *e H . ~ 

ITT T^mVt^l- L' n ~ infl. 

dH^-iuj-f*eE \dH^-iuj-f\x\''^*eE 

Proof: We use the following identities valid for /c-forms in a 3-manifold: 

dF^T) = F*dr], F*{*,rj) = *F*eF*rj, *ceV = c^^^'" *e V- 

Here c is any positive smooth function. If dE = iujii *e H, then 

dE = dF*E = F*dE = F*{iuii *e H) = iuJi*F*e H = iuji *|j|-4e H 
— iuj'jl\x\~'^ *e if. 

The other direction and the other equation are handled analogously. □ 

Next we check how the impedance map transforms. Note that in the form notation, the 
boundary condition in the Maxwell equations corresponds to fixing the tangential 1-form 
*e(i/ A H) on where v — ui dx^ + V2 dx^ + dx^ is the outer unit normal written as a 
1-form. 

Lemma 6.2. Let A he the impedance map in VL with coefficients (/U, 7), and let A he the 
impedance map in Vt with coefficients (/I|5;|~^, 7|5;|~^). If T is a tangential field on dfl, 
then 

A(r) = F*A((F-^)*r). 

Proof: We take p to be a boundary defining function for Q, that is, p is a smooth function 
— > M and Q = {p > 0} , dQ = {p = 0} . Then p = F*p is a boundary defining function 
for Q. The unit normal is the 1-form 

dp F*dp dp ( 2 dp \ \~\2t.* 



—P* 1 — __p* \x\~ — = \x\ F*!/, 

\F*dp\F*{F-^)-e \dp\{F-i)*e \ Mp|e 

using that {F'^Ye — |x|~^e and u — —dp/\dp\e. 
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Let {E, H) be a solution to the Maxwell system in f2 with *e(^ A H)\q^ = T. Then on 
dfl one has 

AE) = \x\^ *e AE) = A E)) = F*(*,(z/ A E)). 

Similarly A H) = F*(*e(z/ A H)), and therefore 

A(f ) = *^{uAE) = F*A{*e{jy A i/)) = F*A{{F-^)*f). 

□ 

We now assume that (fijjjj) are two sets of coefficients satisfying the assumptions of 
Theorem 2, and let Ai and A2 be the corresponding impedance maps. Let Bq be an open 
ball with Q C Bq, and suppose that F C dQ is such that AiT|r = A2T|r for tangential 
fields T supported in T. Assume that Tq = dQ \ T satisfies Fq = dfl fl dBo and Fq 7^ OBq, 
that is, the inaccessible part Fq is part of a sphere. 

We wish to show that fii = fi2 and 71 = 72 in fl. To this end, choose coordinates so 
that Bq = B{xo, 1/2) where xq = (0,0, 1/2), and assume that the origin is not in Q. Let 
Q = F~^{Q) be the image of Q under the Kelvin transform. From Lemma [6.21 we obtain 
that 

Ai(f)|p = A2(f)|p 

for tangential fields supported in F = F~^(F). Since the Kelvin transform maps Fq onto 
a subset of {xs = 1}, we are in a situation where the inaccessible part of the boundary is 
part of a hyperplane. Also the other conditions in Theorem 1 are satisfied; in particular, 
(10. 5p follows from its analog for OBq which reads 

J there exist extensions of 7j and fij into which are 
preserved under the map x F o R o F^^{x), 

where F is the Kelvin transform and R is the reflection (xi,X2,X3) ^-^ (xi,X2,2 — X3). 
Consequently, Theorem 1 implies that 

Jli\x\''^ = Ji2\x\~'^ , 7i|5;|~^ = 72|5;|~^ in f2. 

It follows that fii = ^2 and 71 = 72 in f2, thus proving Theorem 2. 
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